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Abstract. We study the coherent inelastic diffraction of very weakly bound two 
body clusters from a material transmission grating. We show that infernal transitions 
of the clusters can lead to new separate peaks in the diffraction pattern whose angular 
positions determine the excitation energies. Using a quantum mechanical approach to 
few body scattering theory we determine the relative peak intensities for the diffraction 
of the van der Waals dimers (02)2 and H2D2. Based on the results for these realistic 
examples we discuss the possible applications and experimental challenges of this 
coherent inelastic diffraction technique. 

PACS numbers: 36.40.-c, 03.75.Be, 34.50.Ez, 36.90,+f 
1. Introduction 

Since the first observation of Fresnel diffraction of àtoms from a single slit in 1969 
U diffraction of molecular beams transmitted through material devices |J has now 
become subject to intensive studies. New applications involve the use of transmission 
gratings as quantum mechanical mass spectrometers [||, ||, the determination of atòmic 
and molecular elèctric dipole polarizabilities with interferometric precision ||, sensitive 
probing of the interaction of àtoms and molècules with sòlid surfaces || 0] as well as 
fundamental tests of the quantum mechanical nature of complex molècules ||. These 
recent developments all depend on micro-fabricated transmission gratings whose periods 
have now become as small as 100 nm. These small structures allow to scatter collimated 
molecular beams coherently. For helium beams, e.g., far more than 10 diffraction orders 
have been resolved whose intensities extend over more than four orders of magnitude. 

With this great sensitivity diffraction from material transmission gratings can now 
be used for quantitative studies: Several recent experiments have shown that classical 
wave òptics is insufficient to explain the relative magnitudes of the diffraction intensities 
H ||. It turns out that a more general approach based on quantum mechanical scattering 
theory is required to include the van der Waals interaction between the àtoms and the 
grating bars and to account for the finite size of weakly bound clusters in the beams. The 
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highly precise diffraction experiments together with this improved theoretical analysis 
have recently allowed to determine the C3 coemcients of the atom-surface van der Waals 
interaction for a variety of àtoms and molècules ||, as well as the bond length of 
the helium dimer || which characterizes all low energy binary scattering properties of 
helium. 

These transmission grating diffraction studies only involve elàstic scattering. Àtoms 
or molècules scattered from a sòlid surface can undergo infernal transitions. In reflection 
experiments from surface lattices the coherent transfer between H2 molecular rotational 
levels was observed, e.g., in reference | Ï0|| . The possibility of coherent inelastic diffraction 



from transmission gratings was discussed in reference |ÏTJ for the particularly interesting 
example of the excitation of the weakly bound helium trimer 4 He3 whose single excited 
state is believed to be an Efimov state (see, e.g., reference Q)- These studies have 
shown that, quite generally, an inelastic diffraction pattern exhibits peaks that are 
separated from the strong elàstic diffraction màxima with an angular shift depending 
on the energy of the infernal transition and the grating period. The principle of inelastic 
scattering from a material transmission grating, i.e. the (de)excitation of transmitted 
àtoms or molècules, was demonstrated in subsequent experiments using a fine structure 
transition of meta-stable argon àtoms JÏ3| and a vibrational transition of meta-stable 



nitrogen dimers jnj. In these experimental studies, however, the transmission grating 
served only to multiply the intensity through the simultaneous scattering by all the 
parallel surfaces of the bars. The inelastic diffraction peaks implied by the periodicity 
of the grating have not yet been resolved. 

Motivated by the recent progress in inelastic transmission scattering experiments 
in this article we study the coherent (de) excitation of very weakly bound two body 
clusters in the diffraction from a transmission grating. Based on quantum mechanical 
few body scattering theory we discuss to which extent the small momentum transfers 
that the clusters experience when they pass through the slits can excite the low energy 
infernal transitions. We show that the small angular shifts of the inelastic diffraction 
peaks can be well suited to determining the transition energies. The theory is applied to 
the realistic examples of the van der Waals clusters (D 2 )2 and H 2 D 2 . Their level spectra 
are well known both theoretically and experimentally JÏ5, |Tj| [Ï7|, [18]. Furthermore, 



experimental beam sources for these clusters are available. 
2. Diffraction of weakly bound dimers 

In this section we will briefly outline the few-body scattering approach to the diffraction 



of weakly bound two-body systems [[19] that we will apply in this article. For simplicity, 
we denote the two body systems as dimers although, quite generally, the constituents 
may be different àtoms or tightly bound small molècules. These constituents are 
considered as point partides with masses m\ and m 2 that interact with each other 
through a potential V. This potential is assumed to support shallow bound states </> 7 
with negative binding energy Ej. 



Inelastíc Diffractíon and Spectroscopy of Very Weakly Bound Clusters 



3 




Figure 1. Scheme of the material transmission grating and the coordinate system 
used in this article. The bars are equally spaced with a period d. The slit width of the 
grating is denoted by s. The coordinates are chosen in such a way that the z direction 
is parallel to each bar and the grating is periòdic in the y direction. 



The positions of the constituents are denoted by X\ and x 2 . The corresponding 
center of mass and relative coordinates are R = (miXi + m 2 x 2 )/(mi + m 2 ) and 
r = x\ — X2, respectively. The diffracting object for the beam of dimers, i.e. a 
material transmission grating in the applications in this article, is modeled by an external 
potential W(x 1 , x 2 ). For the weakly bound dimers, under consideration, this two body 
potential is given, to an excel·lent approximation, by the sum of the potentials of each 
constituent, W(xi,x 2 ) = W\{ X \) + W 2 (x 2 ). With the diffraction from a transmission 
grating in mind we idealize the diffracting object in such a way that W is invariant with 
respect to translations along the z coordinate (see figure p. This reduces the scattering 
problem, in the center of mass coordinates, effectively to two dimensions. The complete 
Hamiltonian is then given by 

H = + V(r) + W 1 {x 1 ) + W 2 (x 2 ), (1) 

where M = vn\ + m 2 and fi = m 1 m 2 /(m 1 + m 2 ) denote the total and reduced mass, 
respectively. 

A beam of dimers in the state <f)y can be idealized as a stationary energy state 
of the free center of mass motion. For an incident center of mass momentum P' this 
state is given by \P', </y) = \P') where \P') is a plane wave. The total energy 

of a dimer in the beam then consists of the kinetic energy and the binding energy, 
i.e. E' = P' 2 /2M + Ey. Within the usual experimental range of (kinetic) beam 
energies infernal excitations of the constituents should be negligible. As the dimers, 
under consideration, are weakly bound, however, excitations of the relative motion of 
their constituents as well as break-up are possible as long as the center of mass kinetic 
energy exceeds the dissociation threshold. We will study in this article the (de)excitation 
of dimers from the bound state 0y to (f) 1 at beam energies that are large in comparison 
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to the break-up threshold. We will thus assume 

|£ |«(^ 2 + ^ 2 )/2M, (2) 

where E is the ground state binding energy of the dimer. We shall further presuppose 
the diffraction condition, i.e. we assume that the de Broglie wave length of an incoming 
dimer, A' = 2Trh/P', is by far smaller than all length scales set by the diffracting object. 
The characteristic length scales of a transmission grating are the slit width s and the 
width of the bars d — s (see figure 0). Under typical experimental conditions the ratios 
s/X' and (d — s)/X' are of the order of several hundreds or larger. The diffraction 
condition assures that the outgoing intensity is sharply peaked about the incident beam 
direction with a typical width of several mrad. 

The scattering matrix, S, maps the incoming state \P', 0y) onto the outgoing 
state, i.e. the asymptotic state after the diffraction from the potential W. The 5* matrix 
element that describes the (de)excitation of a dimer can be decomposed as 

(P, 7 | S \P', </y> = 5 {3) (P - P')5 lY - 2vri 6(E - E')5(P Z - P' z ) 

x t(P a ,P tf ,0 7 ;J* i* (3) 

where E = P 2 /2M + E 1 is the energy of the outgoing dimer. In reference |Tí| the 
dimer transition amplitude t(P x , P y , 7 ; P x , P y , (fiy) has been determined from a general 
approach to few-body multi-channel scattering theory |2Q| , by means of perturbation 
theory. According to this approach the dimer transition amplitude can be expressed in 
terms of "point particle" transition amplitudes of the constituents: 

tï P (Ptx,Piy; Ap iX} Apiy) = (Ptx,Piy\ Wi \p' ix ,p'iy, +). • (4) 

Here p' ik and pn-, i = 1,2, k = x,y denote the momentum components of a single 
constituent and Ap ik = p ik — p' ik the corresponding momentum transfers. Furthermore, 
P'ixi P'iyi +) • ^ S ^ ne stationary outgoing scattering state that corresponds to the diffraction 
of constituent i with incident momentum {p' ix) p' iy ) from the potential Wi |Ï9| . From 
reference [Ï9[] one then obtains 

t(P x , P y , 7 ; P' x) P' y) </y) = íf p (^P x , ^P y , AP X , AP„) F 77 , (~^P) 



-PP 



H *r ( ^Px, ^P B ; AP X , AP V ) Fw ( ^AP 



2 \M Xí M V1 *' V 77 VM' 
2vriM 



" K t7 (^P,, ^rPy; AP X , ^AP y - , 

P x J \M M M y 

P*> ^ p v-i °> Jï Ap y + %) F W (- ^ AP - °) > ( 5 ) 



2 V M 



where 



1 77 



(p) = J^ r exp(— ip • r/h) 0*(r)0y(r) (6) 

is usually referred to as the form factor of the dimer that corresponds to the transition 
from 0y to 7 . 
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In reference the transition amplitude (^) has been expressed in terms of a 



Fourier transform of a point particle transmission function r. ; 



pp 



*i (Pix,Piy] APix, Ap iy ) = -1 



Pi. 



X 



dï/ exp( 



(2ir) 2 mih 
iAp iy y/h) 



1 



PP 



(pL p\ v \ y) 



1, 2 (see Appendix): 



(7) 



This representat ion is particularly useful in the description of diffraction scattering close 
to the incident beam direction. If, for instance, the interaction of the constituents with 
the grating bars is assumed to be purely repulsive rf p recovers the grating transmission 
function of classical òptics and equation (|7|) becomes the classical Kirchhoff diffraction 
amplitude [22|, chap. 8.5]. 

We shall use equation (|7]) to represent the dimer transition amplitude @ in terms of 
a dimer transmission function: Inserting equation (|7|) into equation (d) and performing 
the momentum integrals shows that, under the diffraction condition and within the 
range of validity of assumption (fj), some terms in equation (|5]) cancel. The dimer 
transition amplitude then assumes the form: 

Per. 

t{Pxi Pyi <Ay! P x i Pyi §i 



— 1- 



X 



(27r) 2 Mh 
dY exp(-iAP y Y/h) 



A _ _dim/ pl pl. y\ 
77 77' \ xi 1 yi 1 ) 



where the dimer transmission function is given by 



t^(P> } P';Y)= /d 3 r 



v' r 



(8) 



(9) 



x r 



pp 



ÜÜP' nUp'-Y 4- ^i\r pp Í^P' Ü^P'-Y-^v 
M M y M J \M M v> M y 



The coordinates R = (X, Y, Z) and r = (x, y, z) in equation 
center of mass and relative coordinates, respectively. 



can be interpreted as 



3. Diffraction from a Transmission Grating 

In this section we will analyze the kinematic diffraction phenomena that are implied by 
the periodicity of a material diffraction grating consisting of equally spaced bars with 
a period d and a slit width s (see figure 0). The grating potential of each constituent 
is assumed to be strongly repulsive inside the bars with an attractive part along their 
surfaces that accounts for the van der Waals interaction of the constituents with the 
material ||. The periodicity of the grating along the y axis reappears in the point 
particle transmission functions rf and thus, by equation (|9]), the dimer transmission 
function T^™ is also periòdic in the center of mass coordinate Y . As a short calculation 
using equation (|9|) shows, this periodicity implies the conservation of the y component 
of the center of mass momentum of a dimer, up to reciprocal lattice vectors[|: 

AP y = n27Th/d, n = 0,±1,±2, .... (10) 
| This conservation law can also be deduced from general symmetry considerations. 
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Furthermore, the conservation of the total energy and the translational invariance of 
the grating along the z axis imply 



P' x l + P?)/2M + R, = (i* + Bj)/2M + P 7 . 



These conservation laws determine, in turn, the angles of the principal màxima of 
the diffraction intensity: The angle of the nth order principal diffraction maximum, 9 n , 
is given, in terms of the momentum P y , by P y = P sin 9 n . In a similar way, P' y = P' sin 9' 
determines the angle of incidence of the beam in the (x, y) plane perpendicular to each 
grating bar. Equations (|Ï0f) and ( |ÏTD then yield 

w \~ l/2 



sin 6 r , 



P' 2 /2M ) 



sin 9' + n 



2nh 

~pül 



(12) 



For elàstic (E y = Ey) diffraction this reproduces the well known formula from wave 
òptics. For inelastic (P 7 ^ Ey) diffraction 9 is shifted to a larger (smaller) angle in the 
case of excitation (de-excitation) and the spacing of the diffraction màxima is increased 



(decreased). We note that equation fll2]) only contains the (de)excitation energy E 1 —Ey 
but no further properties of the bound state wave functions. A physical interpretation 
of equation (O), in terms of refraction of molecular beams, was given in reference fTI 



The break-up of a dimer leads to diffuse scattering angles of the fragments and will not 
be considered further in this article. 

An experimental molecular beam exhibits a finite divergence which leads to a 
broadening of the diffraction màxima. The intensity of a diffraction order is usually 
obtained from the area under the corresponding peak. With this averaging procedure 
the diffraction intensities become independent of the beam properties as long as the 
peaks are resolved. The diffraction intensity of the nth order principal maximum for a 

and éL is then determined, in terms of the S 



transition between the bound states 
matrix (^), through 



J^oc /d 3 P /d 3 P>(P') (P,0 7 |S|P',<M 



(13) 



Here the momentum distribution \ip(P')\ 2 accounts for the beam divergence and C n 
denotes a small cone centered about the nth order principal diffraction maximum. The 
integration over the cone determines the area under the peak. 

To resolve the diffraction orders ip(P') should be sharply peaked about an average 
momentum P' which defines the incident direction as well as the mean velocity of the 
beam of dimers. Under the assumption that the diffraction orders are resolved equations 
(|), © and © yield 

d 2 



oc 



(14) 



where the integral extends only over one period of the grating. The intensities of the 
diffraction peaks are thus determined by the properties of the transmission function 



r. 



dim 



77 
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4. Examples: Diffraction of (D 2 )2 and H 2 D 2 

We will now apply the above results to the van der Waals dimers (D 2 ) 2 and H 2 D 2 . To 
an excel·lent approximation the dimer binding potential V remains unchanged under 
isotopic substitution. The different small reduced masses of (D 2 ) 2 and H 2 D 2 , however, 
lead to pronounced isotope shifts in their bound state spectra. The binding energies of 
both species are of the order of only a few hundred /xeV. Along with (H 2 ) 2 these dimers 
were extensively studied, both theoretically and experiment ally, after finding evidence 
of (H 2 ) 2 in astrophysical planetary spectra. According to their total nuclear spin one 
distinguishes between the ortho (o) and para (p) modifications of H 2 and D 2 . All of 
them may combine to form dimers, leading to a rich infrared fluorescence spectrum (see, 



e.g., references [15 , Ï6| , [Ï7], I8[ , and references [^, for a recent overview of binding 
potentials). 

To facilitate the following discussion we will assume that, by passing the gas 
through a catalytic converter or by a similar experimental technique, the molecular 
beam contains only o-D 2 and p-H 2 . Both molècules exhibit only even rotational 
angular momentum quantum numbers j. Furthermore, we shall assume that the nozzle 
temperature of the beam apparatus is sufficiently low that the majority of the D 2 and 
H 2 molècules are in their rotational ground states j = . 

To determine the dimer transmission function from equation (^) we first have to 
calculate the bound state wave functions 7 (ï") that describe the dimers (D 2 ) 2 and H 2 D 2 . 



To this end we have chosen the semi-empirical potential V(r) by Buck et al |26| which 
is given in analytic form. Its anisotropic contributions will be neglected as we assume 
the constituents to be mostly in spherically symmetric j = states. Then, as usual, 
the stationary two body Schròdinger equation can be expanded in terms of spherical 
harmonies Y/ 11 ^, tp) and one only needs to determine the radial part. A numerical 
integration of the radial equation yields four rotational bound states for (o-D 2 ) 2 and 
three for p-H 2 -o-D 2 , all belonging to the lowest vibrational mode. As higher vibrational 
modes are unbound the radial bound state wave functions Ri(r), as depicted in figure §, 
are determined solely by their end-over-end rotational quantum number l. Within our 



approximations the binding energies of (o-D 2 ) 2 are in agreement with reference JÏ6 . 

The integration over the relative coordinate r = (x, y, z) in the dimer transmission 
function can be performed most conveniently in spherical coordinates (/ d 3 r = 
J oo r 2 dr Jq W d(p Jq sin $d$) . By choosing y = rcosí? the angular integration over ip 
becomes trivial and the substitution a = cos $ yields 



T?™ ml {P' x) P' yl Y) = 5 mm ,V2Ï+ÏVW+ï J o r 2 dr Ri{r)R v {r) J da II™ i 



a) 



pp fm 1 mi , m 2 \ PP fm 2 D/ m 2 ^ m x . 
XTl \M P -M P ^ Y+ M ar ) T > \M P -M^ Y -M ar ^ (15) 



Here the function (et) has been introduced as 



(a) 



1 


(i- 


\m'\)\ (V - 


(m'l)! 


n 




\m'\)\ {V + 


\m'\)\ 



Pl m \a)Pl m \a), (16) 



Inelastíc Diffractíon and Spectroscopy of Very Weakly Bound Clusters 



8 



a 

> 



5^ 



tí 

> 



1=0, E = 

b 


-848 (xeV 


— 1=1, E,= 

b 
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Figure 2. The radial wave functions Ri(r) of (0^2)2 (top) and P-H2-0-D2 (bottom), 
where l denotes the end-over-end rotational quantum number. The corresponding 
binding energies are given in the legend. All bound states belong to the lowest 
vibrational modes. Higher vibrational modes are unbound for both species. The radial 
coordinate is given on a logarithmic scale. The ordinates refer to the binding potential 
V(r) from reference pq| which is shown for comparison. 



where P™ are the associated Legendre polynomials. 11^' (a) is symmetric in the indices 
/ and /'. For identical constituents the symmetry properties of the associated Legendre 
polynomials with respect to reflections at the origin imply a parity conservation selection 
rule for the inelastic diffraction intensity: 

l' + l = even. (17) 

This selection rule can be obtained from equation (^) when one inserts m\ = m<i = M/2 
but may also be derived directly from general symmetry properties of the Hamiltonian 

(0)- 

We assume the dimer bound states in the incident beam to be populated in 
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accordance with the thermal equilibrium weight factors py = (21' + 1) exp(— Ei'/k-^Tb) 
where is the translational beam temperature and k-g is the Boltzmann constant. 
Inserting equation (|Ï5|) into equation (0) and summing over m! we find that the nth 
order diffraction intensity for a transition V— ú is given by 

-J dYexp(-i27rnY/d) Jj 2 àr Ri{r)R v (r) j da Tï% (a) 



/f oc (21 + 1) vv y: 



where, due to the factor 5 mm ' in equation (|Ï5l) , only the orientation quantum numbers 
m' between — mm(l,l') and min(/,Z') contribute to the sum. 

By considering D 2 and H 2 as point partides we have neglected transitions between 
rotational states of the constituents. Due to their small moments of inertia, however, 
the j=0— >2 transition energies are comparatively large, i.e. 22 meV for 0-D2 and 44 



meV for p-H 2 [25]. By keeping the kinetic energy in the beam below these thresholds 
excitation of higher j states can, therefore, be ruled out. 

The diffraction patterns in the following two subsections (figures |] and f|) were 
calculated with experiment ally realizable parameters in mind. The mean beam velocity 
was chosen as v' = P'/M = 500 ms" 1 with a velocity spread of Av'/v' = 8% 0. This 
corresponds to a kinetic energy of 5.2 meV for D 2 (2.6 meV for H 2 ), and a translational 
beam temperature of Tb w 0.4 K in a pure D 2 -beam. We have considered diffraction 
from a typical silicon nitride grating, as characterized experimentally in reference |J, 
with a grating period of d = 100 nm, a slit width of s = 60 nm, and a wedge angle of the 
grating bars of (3 = 5°. The attractive van der Waals interaction between the grating 
and the constituents is accounted for in the point particle transmission functions. Its 
C3 coeflicient (see Appendix) has been determined experimentally for silicon nitride and 
D 2 as C3 = 0.32 meV nm 3 |J. We are not aware of a corresponding experimental value 
for H 2 which we therefore estimated: According to Hoinkes' empirical rule p8| , [|, C3 
is, quite generally, proportional to the static elèctric dipole polarizability of the atom or 
molecule. The polarizabilities of H 2 and D 2 differ by about 1% only. We thus assumed 
the same C3 for both D 2 and H 2 . 

In a typical experiment the diffraction màxima are broadened due to the beam 
divergence and the width of the detector slit. The finite velocity spread of a beam 
leads to an additional broadening which increases with the diffraction order n. These 
mechanisms are all specific to the experimental setup. We, therefore, use an empirical 
model to include the broadening. Previous experimental studies || || have shown that 
the màxima can be represented, to an excel·lent approximation, by Gaussians of the form 
I í , n - \2 / ,,2 



exp i-(e-e n y/wi 



whose widths w n depend on the diffraction order n as 



w n = w Q 



\ 



, , Aw 
1 + ( n 

w 
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Here wq is the angular width of the zeroth order diffraction maximum and Aw accounts 
for the velocity spread. For the beams in the present applications we have chosen 
wq = 3 x 10~ 3 degrees and Aw = 7 x 10~ 4 degrees. Each diffraction maximum is thus 
represented by a Gaussian of this kind with an area I given by equation (|Ï8|). The 
angle of incidence 9' was chosen in such a way that the zeroth order principal diffraction 
màxima of all inelastic transitions do not overlap with the zeroth order elàstic peak. 

4-1- Diffraction of (o-D 2 )2 

In the case of (o-D 2 ) 2 both constituents are identical and the selection rule ([[7]) applies. 
Table shows the allowed transitions and the corresponding excitation energies. Due to 

Table 1. Allowed excitations V— ú for (o-D2)2 and transition energies for the dimer 
binding potential from reference [Era . Also allowed are the reverse processes (de- 
excitation) with negative transition energies. 



l'^l 


E l - E v [ M eV] 


elàstic (/' = /) 







377 


1-V3 


604 



the low translational beam temperature of T& ~ 0.4 K most dimers are in the rotational 
ground state V = 0. The population of the V — 1 state is only of the order of 7%, and 
the initial populations of higher excited states are negligible. 




13.8 13.9 14.0 14.1 14.2 14.3 14.4 
Diffraction Angle 9 [deg] 



Figure 3. Diffraction pattern of a pure (0-02)2 beam including all allowed transitions 
The angle of incidence is 9' — 14°. The intense diffraction màxima are due to 
elàstic scattering. The italic labels indicate the inelastic 0— >2 diffraction orders. Other 
inelastic transitions are masked or too weak to be resolved. The diffraction pattern is 
normalized to the total elàstic zeroth order intensity Iq = Íq° + Iq 1 + Iq 2 + Iq 3 . 
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Figure |3] shows a (o-D 2 )2 diffraction pattern with an angle of incidence of 9' = 14°. 
The intense màxima correspond to elàstic scattering whereas the weak màxima result 
from the inelastic 0— »2 transition. The even weaker 1 — ^-3 màxima are too close to 
the elàstic diffraction peaks to be resolved. According to equation ( |T2"D the spacing of 
the 0—^2 màxima and the shift of the zeroth order with respect to the incident beam 
direction determine the excitation energy E 2 — E = 377/ieV. The hierarchy of the 
inelastic màxima depends on the specific interaction with the grating bars. 



4-2. Diffraction of P-H2-0-D2 

In the case of p-H 2 -o-D 2 the selection rule (|ÏTD does not apply, allowing all transitions 
between the three bound states. Their excitation energies are shown in table |2|. Again, 



Table 2. Allowed excitations V— *l for P-H2-0-D2 and transition energies for the 
dimer binding potential from reference [ p6| . Also allowed are the reverse processes 
(de-excitation) with negative transition energies. 



l'->l 


Ei - E v [^eV] 


elàstic (/' = l) 







170 


0^2 


486 


l->2 


316 



most dimers are in the rotational ground state /' = 0. About 2% of the dimers are in 
the V = 1 state, and the initial population of the V = 2 state is negligible. 




15.8 15.9 16.0 16.1 16.2 16.3 16.4 



Diffraction Angle 9 [deg] 

Figure 4. Diffraction pattern of a pure p-H 2 o-D 2 beam including all allowed 
transitions. The angle of incidence is 9' = 16°. The intense diffraction màxima are 
due to elàstic transitions. The italic labels indicate the resolvable inelastic diffraction 
orders. The diffraction pattern is normalized to the total elàstic zeroth order intensity 

t _ rOO 1 rll 1 7-22 
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Figure [| shows a p-H 2 -o-D 2 diffraction pattern with an angle of incidence of 
6' = 16°. Besides the intense elàstic màxima the inelastic transitions 0—>l and 0— >2 
are resolved. Their spacings determine the transition energies of E\ — Eo — 170/ieV 
and E 2 — E = 486/ieV, respectively. The 0— >1 diffraction màxima are, in general, 
about an order of magnitude more intense than the 0— >2 ones. It is interesting to note, 
however, that the 0— >1 zeroth order is strongly suppressed because the point particle 
transmission functions of H 2 and D 2 are almost equal and almost symmetric. The zeroth 
order vanishes exactly in the hypothetical limit that the van der Waals interaction of 
the constituents with the grating vanishes. Also visible, but weaker because of the small 
initial /' = 1 population, are some màxima of 1— >0 and 1— *2 transitions. 



5. Conclusions 



We have shown in this article how inelastic diffraction from a transmission grating may 
be used to study energy spectra of very weakly bound two body clusters. The transition 
energies are determined by angular shifts of the inelastic diffraction peaks in comparison 
to the strong elàstic peaks. The emergence of these peaks is implied by the periodicity 
of the grating. We have determined their relative intensities from a quantum mechanical 
few body scattering approach for the examples of the van der Waals dimers (D 2 ) 2 and 
H 2 D 2 . 

For these realistic examples the angles of the inelastic diffraction peaks can be 
resolved in present day experiments while their small intensities would require an 
additional gain in sensitivity. This gain may be achieved, for instance, by increasing the 
detection efhciency and, thereby, keeping the background detection rate at a constant 
level. 

Given a suitable detection efhciency the inelastic diffraction technique described in 
this article should be applicable also to those very weakly bound clusters that are not 
easily excited by làser light. A prominent example of this kind, the helium trimer, was 
discussed in reference [11]. Inelastic diffraction from a transmission grating may also be 
used to separate beams of clusters in particular interna! states from the incoming beam. 
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Appendix. The Point Particle Transmission Function 



In reference [^TJ the point particle transmission function was derived in terms of an exact 
scattering solution of the two-dimensional stationary Schròdinger equation as (cf. figure 
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T TP (p' x ,p' y ;y) = expf-vy/AM '^)- 

Here this scattering solution ip(x,y) assumes, at large distances from the grating, the 
asymptotic form of a coherent superposition of the incoming plane wave exp[i(p^x + 
p' y y)/fc\ and an outgoing cylindrical wave. In this appendix we will explain the methods 
we have used to determine T PP (y). We shall focus on off normal incidence of the incoming 
beam (8' ^ 0) as the case of normal incidence on the transmission grating was discussed 

The coordinates as well as the typical trapezoidal shape 
are depicted schematically in figure [Al|. As the grating bars 



in detail in references [[21 , 
of the grating bars 

are assumed to reflect those àtoms that touch the bar walls the wave function (p(x,y) 
vanishes at the surface of the bars. This implies r pp (y) = for the coordinates y at the 
rear of a grating bar. The point particle transmission function T PP (y) in the slits can be 
determined in accordance with the Standard eikonal approximation ||29|| : When, at off 
normal incidence, 6' exceeds the wedge angle f3 a portion t(tan#' — tan/3) of each slit is 
shaded by the adjacent grating bar (see figure |AÏ|). In this case r pp (y) also vanishes in 



the shadow region. For the coordinates y inside the remaining open fraction of a slit the 




Figure Al. The trapezoidal shape of the grating bars with thickness t and wedge 
angle f3, and the direction of the incident beam (dashed line). For 9' > (3 a portion of 
the slit is shadowed (hatched region). 



point particle transmission function assumes the form r pp (y) = exp[i$(y)]. The phase 
shift accounts for the attractive van der Waals interaction between the atom or 
molecule and the grating bars. It is determined by 

*(y) = ~ Jàu W M (u), (A.l) 



where v' = ^p'£ + Py /m is the velocity of the incoming atom or molecule and the 
integration is performed along a straight line in the direction of the beam (dashed line 
in figure |A1| ) . 

The long range van der Waals potential between a pair of àtoms (or molècules) 
separated by a distance R is proportional to — i?~ 6 . To determine the functional form of 
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the attractive interaction between an atom or molecule and the grating one may use a 
pairwise summation of the interactions between the probe particle and the constituents 
of the sòlid. For an idealized flat surface, i.e. when the sòlid covers the half space, 
the pairwise summation implies the well known functional form W att = — C3/L 3 , where 
L denotes the distance between the probe particle and the surface. The coeficient 
C3 depends on the static elèctric dipole polarizability of the atom or molecule and 
the material of the sòlid f| p8| . For D 2 and a silicon nitride grating C3 is known 
experiment ally |J. For H2 C3 has been estimated as explained in Section f|. 

For a grating bar with the geometrical form in figure [Aï| the attractive interaction 
Watt does not assume a simple form. By interchanging the pairwise summation with 
the integral in equation (|A.1|) , however, the phase shift may be determined analytically. 
Numerical studies have shown that only those bars contribute significantly to W a tt that 
confine the actual slit. At off normal incidence with 9' > (3 the phase shift is then given 
by 

C 3 



2/m/Ícos 9 



A4 



en 2 



£f 2 2 



tan 9' + tan (3 



where f n = s/2-y, £ 12 
and £22 = s/2 + y. 



s/2 + í(tan/3 + tan#') 



£,21 ~ ^22 2 

tan 9' — tan/? 
V, 621 = s/2 + í(tan/3 



(A.2) 
tan 9') + y 
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